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Abstract 

A second order finite-difference equation has two linearly independent solutions. It is shown here 
that, like in the continuous case, at most one of the two can be a polynomial solution. The unique¬ 
ness in the classical continuous Heine-Stieltjes theory is shown to hold under broader hypotheses 
than usually presented. A difference between regularity condition and uniqueness is emphasized. 
Consistency of our uniqueness results is also checked against one of the Shapiro problems. An in¬ 
trinsic relation between the Heine-Stieltjes problem and the discrete Bethe Ansatz equations allows 
one to immediately extend the uniqueness result from the former to the latter. The results have 
implications for nondegeneracy of polynomial solutions of physical models. 
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I. INTRODUCTION 


Let A{x) and B{x) be given polynomials of degrees m + 1 and m, respectively. The 
snbject of the classical Heine-Stieltjes theory is to determine a polynomial V(x) of degree 
m — 1 snch that the second-order differential equation 


A(x)y'' -h 2B(x)y' + V(x)y = 0 


( 1 ) 


las a solution which is a polynomial of a preassigned degree n jl-5|. Assume with Stieltjes 
that A(x) has real unequal roots, 


A(x) = (x — ao)(x — ai)... (x — am), ag < ai < ... < am, 


( 2 ) 


and 


B(x) 


Po Pi 


+ ... 


Pv 


> 0, z/ = 0,1, 2,..., m. 


( 3 ) 


A{x) X — ag X — ai x — a._^ 

This is equivalent to the assumption that the zeros of A{x) alternate with those of B{x) and 
that the leading order coefficients of A{x) and B{x) have the same sign. Under the above 
conditions, the basic properties of polynomial solutions are 0.0,00: 


there are exactly 


^ n.m. 


n + m — 1 


n 


polynomials V{x), which are called van U/ecA: polynomials 


001 . 


Eq. ([T]) cannot have two polynomial solutions linearly independent of each other. 


If 1 / is a polynomial solution, y ^ 0, then y 0 a.t x = a,^. 


All the zeros of y are distinct. 


The zeros of y lie in the interval [oq, On 


The case with m = 1 corresponds to the hypergeometric differential equation, while the case 
with m = 2 corresponds to the Heun equation j^. For m < 3 polynomial solutions mostly 
characterize QES models lCHl6l|. although not all such polynomial solutions are exhausted 


by the QES models 


14| . General (extended) Heine-Stieltjes polynomials were often studied 
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in connection with a special Lipkin-Meshkov-Click model corresponding to the standard 
two-site Bose-Hubbard model 0 , 0 . 

Compared to the continuous case of Eq. ([T]), very little is known about general properties 
of polynomial solutions of a linear homogeneous second-order hnite-difference equation 


g{x)A‘^y{x) r{x)Ay{x) -\- u{x)y{x h) = 0, 


where the hrst difference quotient of y{x), or Norland’s operator A 

h 

in usual sense 


HQ, 


(4) 

is dehned here 


Ay{x) = 


y{x + h) -y{x) 


The hnite-difference equation (jl]) can be disguised in further equivalent forms 

g{x)A^y{x) -f [r(x) -|- hu{x)]Ay{x) -f u{x)y{x) = 

h h 

g{x)y{x 2h) + [hr{x) h?u{x) — 2g{x)]y{x -\- h) [g{x) — hr{x)]y{x) = 0. (5) 

(The hrst one follows on making use of the identity ay{x -\- h) = haAy[x) -|- ay{x).) Last 

h 

but not the least, if y{x) = n^=i(^ ~^j) ^ polynomial solution, then (|5]) leads at any 

zero Xk of y{x) to a discrete Bethe yfnsafe equation (cf. Sec. 5 of Ref. [19|) 


]Tj=ii.^k - Xj + h) ^ hr{xk - h) - g{xk - h) 


- h) 


g{xk - h) 


( 6 ) 


The motivation to study polynomial solutions of hnite-diherence equations has got a boost 
after it was demonstrated that physical models with a discrete nondegenerate spectrum can 
be characterized in terms of orthogonal polynomials of a discrete variable and their weight 


function 


20 


25| . The latter applies to all problems where Hamiltonian operator is a self- 


adjoint extension of a tridiagonal Jacobi matrix of dehciency index (1,1) 2^. For instance 
a displaced harmonic oscillator can be characterized in terms of the classical Charlier poly¬ 
nomials and the Rabi model by a norm preserving deformation of the Charlier polynomials 
221 . l23| . Some earlier applications of classical discrete polynomials in physics not related 
to Lanczos-Haydock scheme 20, l2l|] have been given by Lorente 27(]. He showed that the 
respective orthonormal Kravchuk and Meixner functions are related to a quantum harmonic 
oscillator and the hydrogen atom of discrete variable, and that the Hahn polynomials are 
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related to Calogero-Sutherland model on the lattice. 


Unfortunately, only the hypergeometric case m = 1, where the polynomial coefficients 
g{x),r{x),u{x) have degrees 2,1,0, respectively, has been studied exhaustively within the 
realm of c/assfca/orthogonal polynomials of a discrete variable 28l-l3 1L Generalized Bochner 


theorem for hnite-difference equations has been dealt with in Ref. 
forward has been achieved by Turbiner 
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32[ |. An important step 


13l | within the realm of quasi-exactly-solvable 


(QES) equations jl0l4l5l|. The latter yield a specihc subclass of hnite-difference equations 


where the polynomial coefficients g{x),r{x),u{x) have degree at most four. 


The motivation of present work is to translate the properties of the classical continuous 
Heine-Stieltjes theory into the realm of hnite-difference equations. As in the continuum 
case, a second order hnite-diherence equation (jl]) has two linearly independent solutions 
for a hxed triplet of polynomial coefficients g{x),r{x),u{x). Here we derive the conditions 
under which two linearly independent polynomial solutions of Eq. (jl]) are forbidden, i.e. 
the polynomial solutions of general second-order hnite-diherence equation (|1|) are unique 
(cf. Theorems 1 and 2). As a by product, an h-analogue of AbeFs theorem for the Heine- 
Stieltjes problem is derived, which yields an explicit analytic expression of hnite diherence 
Wronskian, or Casoratian, Wh{x) in terms of a rational function involving products of gen¬ 
eralized gamma function r/i(a;) in both numerator and denominator. A comparison with 
the classical hypergeometric equation is provided in Sec. Ill Al Using an intrinsic relation 
between the Heine-Stieltjes problem problem (|1|) and the discrete Bethe Ansatz equations 
(l6ll . the uniqueness result is extended from the former to the latter in Sec. IllB1 The results 
are discussed from diherent perspectives in Sec. IIIIl Sec. IHIAI shows that the uniqueness in 
our sense ensures uniqueness even if the regularity condition for the Hahn class of hypergeo¬ 
metric orthogonal polynomials (cf. Sec. 2.3 of Ref. 3l|]) does not preclude two polynomial 
solutions. A comparison with one of the Shapiro problems is discussed in Sec. IHIBl There 
are two basic ways how to make use of the uniqueness theorems for Heine-Stieltjes polyno¬ 
mials. First, they yield a straightforward proof of the nondeqeneracy of QES levels which 


yield the so-called exceptional spectrum of physical models 
which is more involved by other means (cf. Refs. Q, 
theorems in certain exceptional cases - cf. Sec. IIIIl 
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RHQ). 


13 


, lisl, l2]J, I25I , the proof of 


Second, they serve as no-go 
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II. UNIQUENESS 


For each a; G M one can define the lattice Ah{x) := {x + kh \ k G Z}. For a given 
Xo G M the second order finite-difference equation (E]) is seen to connect the values of y{x) 
at the points of A/j(xo). The function constant on each Ah{xo) is called h-periodic function. 
Two functions yi and y 2 are called linearly dependent in a finite-difference sense if there 
are h-periodic functions Ci and C 2 such that Ci{x)yi{x) -|- C 2 {x)y 2 {x) = 0. Otherwise the 
functions yi and y 2 are called linearly independent. We shall use repeatedly the following 
elementary argument; If y{x) is known to be a polynomial of degree not larger than N and, 
at the same time, to vanish in at least -|- 1 different points, then y{x) = 0. In what follows 
we shall consider the second order finite-difference equation (E]) with polynomial coefficients 
only. Using the argument, one finds immediately that: 

• (PI) If a polynomial y{x) solves equation ([5]) on an infinite subset of A/i(xo), then 
y{x) solves it for all xq G M. 

• (P2) If a linear combination Ciyfix) 6 * 22 / 2 ( 2 ^) of two polynomials vanishes on an 
infinite subset of A/i(xo), then it vanishes for all x G M. 


The latter implies that linear dependence of two polynomial solutions yi and 1/2 in a finite- 
difference sense reduces to the linear dependence in conventional sense, i.e. with Ci and C 2 
being independent of x. 

The following theorem, and Theorem 2 below, encompass all quasi-exactly-solvable equa¬ 


tions on a uniform 
discrete variable 


28 


inear-type lattice 10l-ll3| and all classical orthogonal polynomials of a 




Theorem 1: Let the second-order finite-difference equation (j4]) has polynomial coeffi¬ 
cients such that g{x) and g{x) — hr{x) have real roots. 


9{.x) 

= (x - bo){x 

- 61 )... (x - b^fi, 

V 

V 

0 

.0 


g{x) — hr{x) 

= (x — ao)(x 

fll) . . . (x Um)) 

(Zq . 

• • 


(7) 


For each root aj define a uniform lattice A^^. := {aj kh\k G No}, which extends to the 
right of the root aj. It is not excluded that Oj^, bi G A^^.^ for bi > aj^. Assume further 
that there is at least a single A^^. which does not contain any root of g{x). Then Eq. (j4]) 
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cannot have two polynomial solntions yi and 1/2 linearly independent of each other. 


Proof. 


For any two fnnctions m and 1/2 the Leibniz's theorem of hnite-difference calcnlns (pp. 


34-35 of Milne-Thomson 


18l |l implies 


A[yi{x h) ■ Ay2{x) - Ayi{x) ■ y2ix + h)] = y^{x + h) ■ A^y2{x) - A^yi{x) ■y2{x + h). 

h h h h h 

Hence for two nontrivial solntions yi and y 2 of the hnite-difference eqnation (jl]) we have 
g{x)A[yi{x + h) ■ Ay2{x) - Ayi{x) ■ y2{x + h)] 

h h h 

+r{x)[yi{x + h) ■ Ay2{x) - Ayi{x) ■ y2{x P h)] = 0 . (8) 

h h 

The latter is of the form 


AX(a;) =-^X(x) or X{x + h) = X{x) = R{x)X{x), (9) 

h g[x) g[x) 

where X stands for the sqnare bracket in Eq. ([ 8 ]), which can be identihed with a hnite 
difference Wronskian, or Casoratian, 33 1 


Wh{yi,y2}{x) : = 


yi{x + h) y2{x + h) 

1 

yi{x) y2{x) 

Ayi{x) Ay2{x) 

h h 

“ h 

yi{x + h) y2{x + h) 


= yi{x) ■ Ay2{x) - Ayi{x) ■ y2{x). 

h h 


( 10 ) 


The hypotheses of Theorem 1 determine R{x) as a rational fnnction with zeros and poles 
on the real axis 

g{x) - hr{x) _ ]\%o{.x - aj) 


R{x) := 


( 11 ) 


nr=o(^ - bi) 

Now if Aaj does not contain any zero of g{x), i.e. R{x) is not singnlar on , then, in 
virtne of R{aj) = 0, the hrst-order recnrrence dH]) implies X{x) =0 for all x G Aa^+h- In 
other words, for each Xg G Aa.+h there are C'i(xs) and C 2 {xs) not both zero, snch that 
Ci{xs)yi{xs) R C 2 {xs)y 2 {xg) = Ci{xs)yi{xg + h)+C 2 {xs)y 2 {xs + h) = 0. Taking = Oj + h, 
the linear combination y{x) := Ci{xs)yi{x) -|- C2{xs)y2{x) is a solntion of Eq. (jl]) on Aa^+h 
which satishes y{xa) = y{xs -|- h) = 0. Considering the latter as the initial valnes of the 
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Cauchy problem for the recursive form ([5]) of Eq. (jl]), one has y{x) = 0 on Because 

g{x) 7 ^ 0, the solutions of the Cauchy problem for Eq. ([5]) are uniquely determined by the 


initial values 


33| . and hence Ci and C 2 are constants on entire A^a-+h- In virtue of the 


elementary argument (P 2 ), the linear combination Ciyi{x) + C2y2{x) vanishes for all x G 
i.e. y\ and y2 are linearly dependent in the conventional sense. 


Remark: On considering equation ([9]) as a downward recurrence X{x) = R~^{x)X{x+h), 
an alternative version of Theorem 1 follows which guarantees the uniqueness, provided that 
there is at least a single which does not contain any root of g{x) — hr{x). Here 
is dehned for each root bi as a uniform lattice which extends to the left of the root bi, 
Abi := {bi - kh\k e No}. 

Theorem 2 : Let us consider the second-order hnite-difference equation (jll) with the 
polynomial coefficients as in Theorem 1 . Assume further that there is at least a single Aa- 
which contains more roots (e.g. the single root aj) of g{x) — hr{x) than the roots of g{x) 
[e.g. none of the roots bi of g{x)]. Then Eq. (jl]) cannot have two polynomial solutions yi 
and j /2 linearly independent of each other. 


Before giving the proof of Theorem 2, it is expedient to provide an h-analogue of AbeFs 
theorem which yields an explicit analytic expression of X[x) in terms of a rational function 
involving products of T/i in both its numerator and denominator. The h-extension of the 
gamma function T^fx) is introduced through the functional equation Th{x + h) = xVh{x) 
(cf. sec. 9.66 of Ref. 181: Appendix iBll . 


Lemma 1 : For any rational R{x) of the form flTT]) . the solution X{x) of the hrst-order 
hnite-difference equation ([9]) is either identically zero or 


n^n ~ CLj) 

X (x) = const X —- 

Eon(a:-6,) 


( 12 ) 


Provided that the ratio k of the leading polynomial coefficient of g{x) — hr{x) to that of 
g{x) is K 7 ^ 1, the r.h.s. of Eq. (IT^ will acquire an additional multiplication factor [cf. Eq. 
and becomes 

,UT=o'^h{x - aj) 


X{x) = const X fC-Cm, 




(13) 
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Proof. 

First, Eq. (|9]) is recast as 


h h \ g[x) 


which has the form of the hrst-order hnite-difference equation flAip . Its solution can be 
expressed in terms of Norlund’s principal solution [^, 18|, an elegant, but nowadays largely 
forgotten, tool of integrating finite-difference equations (see Appendix [A] for a brief summary 
and dehnition), as 

'g{t) - hr{t)\ 


X (x) = exp 


^ 1 


In 


gif) 


-} 


At 

h 


(14) 


Note in passing that use of a partial fraction decomposition ([3]) of the fraction in the 
integrand in the exponent of flT^ . as in the continuous case of Stieltjes and further 
elaborated in Sec. 6.81 of Ref. would not bring us any further. Instead it is expedient to 
substitute the respective products (El) into Eq. flT^ and use the logarithm there to split the 
resulting ratio into a sum of individual logarithms ln{t — aj) and — ln{t — bi) corresponding to 
the roots in Eq. (E]). Each such a logarithm term integrates to a corresponding generalized 
gamma function Th (cf. Eq. (IB2p of Appendix [Bl sec. 9.66 of Ref. |l8|). The latter recipe 
enables one to express (IT) as in Eq. flT^ . The transition from (ITT)) to f[T3|l is similar to 
that used by Lancaster 28|] in arriving from his Eq. (29) to his Eqs. (30-33). 


Proof of Theorem 2: 

If X{x) of two linearly independent solutions in Eq. ([9]) is not identically zero, the 
hypotheses of Theorem 2 imply that X is necessarily singular for some its argument value, 
which is impossible if yi and y 2 are polynomials. Indeed, the hypotheses of Theorem 2 
ensure that there is at least a single Aa^ which contains more roots (e. g. the single root 
Qj) of g{x) — hr[x) than the roots of g{x) (e. g. none of the roots bi of g{x)). Unless X{x) 
is identically zero. Lemma 1 determines the analytic form of X{x) to be either (IT2|) or (IT3j) . 
Now Th{x — aj) has a simple pole at a; = aj (cf. Appendix iBl) . If there is Oj < Ok E A^^., then 
also Th{x — Ok) has a simple pole at a: = Oj. If there is bi G Aa-, some of the simple poles 
of Th{x — Oj) and Th{x — a^) at a; = Oj in the numerator on the r.h.s. of Eq. flT3l) could 


8 














be canceled by the simple pole of the rft(a; — bi) at x = aj in the denominator on the r.h.s. 
of Eq. flT^ . Nevertheless, the hypotheses of Theorem 2 gnarantee that at least one of the 
simple poles of T/j’s in the nnmerator is not compensated by the simple pole of r/j(x — bi) in 
the denominator. Then X{x) tends to inhnity for x ^ aj. However, as a discrete Wronskian 
of two polynomial solutions, X{x) cannot tend to infinity at any hnite x G M. Of course, the 
latter does not hold for general nonpolynomial solutions. Thus, as in the continuum case of 
Sec. 6.81 of Ref. we have a contradiction, unless, of course, X = 0. 


A. Classical hypergeometric equation 


As an example, consider the classical hypergeometric equation 28l-l31| 


(ax^ + 6x + c)X^y{x) + {dx + f)Ay{x) + Xy{x + h) = 0. (15) 

h h 

A necessary and sufficient condition for the existence of a polynomial solution of Eq. (IT^ 
is that a characteristic polynomial. 


9{z) := az{z — 1) + dz + X, 


has a non-negative integer root (cf. the n = 2 case of Theorem 2 of Ref. 28(]). If there is a 
polynomial solution of degree n, then 6{n) = 0. The latter is equivalent to 


A -|- nd -|- n{n — l)a = 0, or A„ = —n(n — l)a — nd, n = 0,1,2,.... 


( 16 ) 


Eq. (IT^ is a special case of the eigenvalue problems for the Hahn class of orthogonal 


polynomials 


m, 


3l[ |. In the latter case the regularity condition says that all eigenspaces of 


the hypergeometric eigenvalue problem are one dimensional if and only if A„ 7^ A/ for I 7^ n 
in the set of numbers {A„}))Tq dehned by Eq. flT^ . or if and only if a[n] -|- d 7^ 0 (cf. Sec. 
2.3 of Ref. M])- Here [—1] = —1/q, [0] = 0, [n] = n > 1, and q G M\{ —1,0} is 

the Hahn parameter (for the uniform linear lattice in our case q = 1 and [n] = n). However, 
the regularity condition does not exclude the corresponding eigenspace to be, for instance. 
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two dimensional for = Xi with I 7^ n. The latter is precluded by the following Corollary. 


Corollary: Polynomial solutions of the second-order hnite-difference hypergeometric 
equation Eq. flT^ are nondegenerate, i.e., for a given eigenvalue A there is at most a single 
solution to Eq. flT^ . 

If d = —ka, then Xn = —n{n — k — l)a and Xn may equal A; for some I 7^ n. For instance, 
Xi = Xk = k for n = l,k. \i k > 1 there is thus, under the hypotheses of Theorems 1 and 2, 
no polynomial solution of degree k. 


B. Discrete Bethe Ansatz equations 

Using an intrinsic relation between the Heine-Stieltjes theory and the discrete Bethe 
Ansatz equations one can immediately arrive at the following result. 

Theorem 3 : Provided that the Heine-Stieltjes problem has unique polynomial solution, 
the corresponding discrete Bethe Ansatz equations (El) have also a unique polynomial solu¬ 
tion up to permutations of zeros x^s. 

Proof: 

A solution y{x) = YTj=oyj^^ to the discrete Bethe Ansatz equations ([6]) implies that the 
second-order difference equation (jl]) is satished at the n points xi,X 2 , ■ ■ ■ ,Xn- The necessary 
condition that y{x) = solves Eq. (jl]) is the vanishing of the leading nth degree. 

The latter requires that the sum of the coefficients of the leading degree of the polynomials 
g{x), [hr{x) -|- hfu^x) — 2g[x)], and [g{x) — hr[x)] in the recurrence form ([1]) of Eq. (jl]) 
vanishes. In the hypergeometric case this is the condition flTbl) . If the polynomial coefficients 
of Eq. (jlD are assumed to satisfy the necessary condition, the l.h.s. of Eq. (jl]) becomes a 
polynomial of one less, i.e. (n — l)th, degree. By the elementary argument (e.g. leading to 
PI), if a polynomial in x of degree u — 1, that can vanish only at n — 1 different points, 
vanishes at the n distinct points Xi,X 2 , ■ ■ ■ ,x„, then it must vanish identically. Thus the 
l.h.s. of Eq. (jlD vanishes identically. This leads to a second-order difference equation whose 
polynomial solutions are unique. 
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Ismail et al have earlier shown that the solution to the discrete Bethe Ansatz equations 
(1^ with the right-hand side derived from the Meixner Mn{x] /3, c) and the Hahn polynomials 
Qn{x-, a, / 3 , N) are unique up to permutations (^. Sec. 5 of Ref. 1^). Theorem 3 extends 
the results of Ismail et al (cf. Sec. 5 of Ref. [l9|) to the general case. Some special cases 
when the uniqueness may break down are discussed in Sec. IIIIBI 


III. DISCUSSION 


Our uniqueness theorems encompass all quasi-exactly-solvable equations on a uniform 


linear-type lattice 


ljl3| 


and all classical orthogonal polynomials of a discrete variable 


28 


31| . The hypotheses of our uniqueness theorems look rather different from those in the 
classical continuous Heine-Stieltjes theory [l-5|. In the hnite-difference case, the respective 
g{x) and g{x) — hr[x) can be identihed as the coefficients of y{x + 2h) and y{x) in the 


recurrence form (|5]) of Eq. (jlj). Unlike the continuous case of Refs. [2|, l5| (i.e. with A in Eq. 

h 

dl]) replaced with ordinary derivatives [cf. Eq. ([I])]), one does not assume that the zeros of 
g{x) alternate with those of r{x) [cf. Eqs. (|2]), ([3])]. The hypotheses of Theorems 1 and 2 
are also silent about relative degrees of the polynomial coefficients g{x),r{x),u{x). 

However, the above differences are mostly only apparent, until one realizes that already 
in the classical continuous Heine-Stieltjes theory the assumptions that (i) the zeros of A{x) 
alternate with those of B{x) and that (ii) the leading order coefficients of A{x) and B{x) 
have the same sign, are not necessary for the uniqueness of solutions. Indeed, one can 
multiply both the numerator and denominator in B{x)/A{x) on the l.h.s. of Eq. ([3]) with 
the same polynomial factor (x — 7/)”*, > 1, without changing the r.h.s. of Eq. (|3]), and 

hence the reasoning leading to the uniqueness. It is also not necessary that all > 0 as in 
Eq. ([3]). (The latter has been recognized as late as 2000 by Dimitrov and Van Assche 34|.I 

A broader sufficient condition for the Wronskian W{yi,y 2 } to diverge to inhnity is that 
there is merely at least one u such that p,y > 0 and is different from all other 6^’s. The 
latter points could be illustrated for a continuous hypergeometric analogue of Eq. flTHl) . 


(ax^ -\-bx -\- c)y"{x) -|- {dx f)y'{x) -\- Xy{x) = 0. 


(17) 
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A. Regularity condition vs nniqneness 


The regularity condition of the eigenvalue problems for the Hahn class of orthogonal 
polynomials does not answer what happen if A„ = A/ for I ^ n in the set of numbers 
{A„}^g dehned by Eq. ffTSj) . Will the eigenspace corresponding to A„ = A; be zero-, one-, 
or two-dimensional? The question of uniqueness and existence of the polynomial solutions 
of the hyperqeometric equation flTTl) reduces to solving Lesky’s downward TTRR (cf. Eq. 
(3) in Ref. 351) 


{n - k)[{n + k - l)a d]ank = {k + l)[{k + 2 )can,k +2 + {kb + f)an,k+i] 


(18) 


for the coefficients Qnk of the polynomial solution of the nth degree, 


2/(x) Qjfin^ dfi^n—l^ “h • • • H“ dnO 


(19) 


The TTRR runs downward for /c = n — 1, n — 2,..., 0, with the initial condition an,n+i = 0. 
Without any loss of generality one can assume Onn = 1- With the initial conditions on 
o-n,n+i and Onn being hxed, any other not linearly dependent solution has to have an,n+i 7^ 0 
for Wh{x) [see Eq. (ITOD ] of the Cauchy problem for the TTRR flTSD to be nonzero. This 
is impossible for a polynomial solution of the nth degree, which implies uniqueness of the 
polynomial solution of the nth degree, provided it exists (i.e. {n + k — l)a + d 0). 

The condition ffT6|) is valid both in the continuous and discrete cases. Thus for d = —ka 
some of Xn may equal Xi also in the continuous case (e.g. Ai = Aa, = for n = l,k). 
Let [a:] denotes the smallest integer not less than x, or the ceiling function. Then, unless 
some additional conditions are satished, Lesky’s TTRR flTSll does not have any solution for 
d = —ka and n G [[(fc + 2)/2],/c -|- 1]. Obviously the uniqueness of polynomial solutions 
persists even though the above assumption (ii) is not satished. 


B. Shapiro problem 

The additional conditions under which Lesky’s TTRR flTSl) has a solution for any degree 
n and when the uniqueness of polynomial solutions breaks down are formulated separately 
for k even and odd. The latter is related to the problem of describing when a linear ordinary 
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differential equation with polynomial coefficients admits at least 2 polynomial solutions, 
which is the hrst of hve open problems listed by Shapiro jb]. An exhaustive answer in the 
special case 2B{x) = —A'{x) has been obtained by Eremenko and Gabrielov 36|. The 
following discussion is limited to the hypergeometric equation f[T7|) but is not constraint to 
2B{x) = —A'{x). 

For d = —ka and even k = 2t > 0 (Lesky’s special case 2), uniqueness persists unless 
/ = —tb. Then dx + f = —2tax — tb, or 2B{x) := dx + f = —tA\x) in the notation 
of Eq. ([T]), and hence all the residues pj of the ratio B{x)/A{x) = —tA'{x)/[2A{x)] in 
Eq. (|3]) are necessarily negative. For odd k = 2t — 1 > (Lesky’s special case 3), the 
ratio B{x)/A{x) = —t/{x — oq) — (t — l)/{x — oi), i.e. none of the residues pj of the ratio 
B{x)/A{x) is positive. Thus not just any algebraic dependence of A{x) and B{x) but only 
a particular one jb] leads to that the uniqueness of polynomial solutions ceases to hold and 
there are possible two linearly independent solutions of the continuous Eq. flTTl) for the same 
value of A. 


IV. CONCLUSIONS 


We have established sufficient conditions (Theorems 1 and 2) for the uniqueness of po¬ 
lynomial solutions of second order hnite-difference equations. They encompass all classical 


orthogonal polynomials of a discrete variable 
on a uniform linear-type lattice 
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28l-l3l| and all quasi-exactly-solvable equations 


13| . An h-analogue of Abel's theorem for the Heine- 


Stieltjes problem was derived, which yields an explicit analytic expression of hnite difference 
Wronskian, or Casoratian, Wh{x) in terms of a rational function involving products of gener¬ 
alized gamma function Th{x) in both numerator and denominator. The latter was facilitated 
by Norland’s principal solution'^ [l^, 18|. It suffices to know Norlund’s principal solution2] 
only for a constant [cf. Eq. flAbp ] and a logarithm [cf. Eq. flB2p ] to deal with a large set 
of hnite difference problems (e.g. Ref. j^). Using an intrinsic relation between the Heine- 
Stieltjes problem (|1]) and the discrete Bethe Ansatz equations ([b]). Theorem 3 extended the 
uniqueness of polynomial solutions of the discrete Bethe Ansatz equations of Ismail et al 
(cf. Sec. 5 of Ref. [^) to the general case. The uniqueness in the classical continuous 
deine-Stieltjes theory was shown to hold under broader hypotheses than usually presented 
, y, 1^. A difference between the regularity condition and uniqueness was emphasized. 


13 



















An extension of the results to a general lattice and a second-order finite-difference equation 


with A being replaced by the more general Hahn operator 

h I I 


29 


3l| is dealt with 


in a 


forthcoming publication 37|]. An open question remains if it is possible to translate also 
the remaining properties of the classical continuous Heine-Stieltjes theory into the realm of 
hnite-difference equations. 
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Appendix A: Norlund’s principal solntion 


That particular solutions of the given equation 

Au{x) = 0(x), (Al) 

h 

always exist is seen (in the case of the real variable) by considering that u{x) being arbitrarily 
dehned at every point of the interval 0 < x < h, the equation dehnes u{x) for every point 
exterior to this interval. The expression 

f{x) = A — h[(j){x) + (j){x + h) + (j){x + 2h) + (f)[x + 3h) + ...] 

OO 

= A — h (f){x + sh), 

s=0 

where A is constant, is a formal solution of the difference equation, since 


f{x + h) = A — /i[0(x + h) + 0(x + 2h) + 4>{x + 3h) + ...], 


and therefore f{x + h) — f{x) = h(j){x). However, such solutions are in general not analytic. 
Norland 17| has succeeded in defining a principal solution which has specially simple and 


definite properties. In particular, when (j){x) is a polynomial so is the principal solution. If 
for A we write (j)(t)dt, and if this infinite integral and the infinite series both converge. 
Norland dehnes the principal solution of the difference equation, or sum of the function (j){x), 
as 

£ POO 

(j){z) Az = / 0(t) dt — h ^ 0(x + sh). (A2) 


s=0 


The principal solution thus dehned depends on an arbitrary constant c. As an example. 


consider 


17|,Q 


Au{x) = e 


X and h being real and positive. Here 


F{x) = 


/ OO 

dt — h 

s=0 


—x—sh 


= e-"- 


he 


1 — e 


-h' 


(A3) 
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after evaluating the integral, and summing the geometrical progression. 


The necessary and sufficient conditions for the existence of the sum F{x) as dehned above 
are the convergence of the integral and of the series. In general, neither of these conditions 
is satished and the dehnition fails. In order to extend the dehnition of the sum, Norlund 


adopts an ingenious and powerful recipe. This consists in a regularization of 0(x 
parameter fi (> 0), say 4)(x,u), which is so chosen that (see Chapter III of Ref. 


with a 


also Chapter VIII of Ref. 


17|; see 


13 ) 


• (i) ^ 0(^); 

• (ii) r ' (f){t) dt and both converge. 

For this function 0(x,/r), the difference equation 

Au{x) = (A4) 

h 

has a principal solution, given by the dehnition flA2p . 

pOO ^ 

F(x,/u)= / (pit, ^)dt — h 0(a; + sh, /i). 

s=0 

If in this relation we let /r —?• 0, the difference equation flA4H becomes the difference equation 
flAip and the principal solution of the latter is dehned by 

Fix) = lim Fix, /i), 

fi^O 

provided that this limit exists uniformly and, subject to conditions (i) and (ii), is independent 
of the particular choice of 0(a:,/x). When the limit exists 0(x) is said to be summable. 

The success of the method of dehnition just described depends on the diherence of the 
inhnite integral and the inhnite series having a limit when /i —)■ 0. Each separately may 
diverge when /x = 0 and the choice of 0(x, /x) has to be so made that when we take the 
diherence of the integral and the series the divergent part disappears. It has been shown 
that, for a wide class of summation methods, the result is independent oi the method adopted. 
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A convenient practical choice is 





= lim 

/i—^0 


dt-hJ2^i^ + ^ , 


(AS) 


s=0 


where p > 1, q > 0, such that for A(a;) = x^(hix)'^ this limit exists. Norlund’s recipe 
can be seen as a two-parameter extension of the single-parameter Lindelof and Mittag-Leffler 
methods of summing divergent series 


The latter belongs to the so-called analytic and 


regular summahility methods [38|,l39(]. If applied to a power series (i) it yields the value equal 


to that obtained by an analytic continuation of the series beyond the radius of convergence 
anytime the limit exists, (ii) pro vided that the sum converges for p = 0, the limit p —)■ 0 


yields the very same sum 


38 


39|. 


As a simple illustration, consider 


Au{x) = a, 

h 

where a is constant. The series a -|- a -|- a -|- ... obviously diverges, but for /i > 0 


OO 


/ ae dt, E ae 

s=0 


—fi(x+sh) 


both converge if h is a positive real number, so that we can take A(x) = x, i.e. p = 1, g = 0. 
Hence 



ae dt — h ae 


ae ahe 


/i 


lim 

/i—^0 


lim 

fi —^0 


/i—^0 


lim 

f-L —^0 


h 

a\X -c - - 


—^(x+sh) I 

s=0 J 


1 - 


ae-^" 

ae~^^ 

-f _ e-M/i _ 


/r(l — 

ph — — ph- 

- 

to 

1 

1 

1_1 


[ph + ...] 



(A6) 
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which is the principal solution. It should be noted that both the integral and the series 
diverge when /r = 0 . 


Appendix B: The generalized Gamma function 

Following sec. 9.66 of Ref. [^, if we define the function rft(x) by the relation 

h.lnr/i(a;) = ^^InzA^ + hln \/2ti/ h, 


we have by differencing 


and hence 


hA\nTh{x) := In ^ = Inx, 

h Th{x) 


Th{x + h) = xTh{x). 


(Bl) 


(B2) 


(B3) 


Thus, if n be a positive integer, Th{nh + h) = K^nWhih). r/j(a;) can be related to the 
conventional r(x) through 


lnrft(a:) = lnr(x/h) + — (x — h) Inh, 


or 


Th{x) = T{x/h) exp ( ^ . In h 1 . 


Using the above relation one finds Th{h) = 1, and for any positive integer n > 0 


Th{nh + h) = h^nl. 


The formula (sec. 9.66 of Ref. 



1 


7-ln h „ 

= e h 


n 

S=1 


X 

sh 


(B4) 


where 7 ~ 0.5772 is the Euler-Mascheroni constant, shows that l/r/i(x) is an integral 
transcendent function, with simple zeros at the points 0 , —h, — 2 h., —3h,..., and therefore 
that r/i(a;) is a meromorphic function of x with simple poles at the same points. 
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